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ABSTRACT 

Purpose: Most applications in the field of medical image processing require precise 
estimation. To improve the accuracy of segmentation, this study aimed to propose a novel 
segmentation method for coronary arteries using computed tomography angiography (CTA) 
data sets to allow for the automatic and accurate detection of coronary pathologies. 
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Methods: The proposed segmentation method included 2 parts. First, 3D region growing was 
applied to give the initial segmentation of coronary arteries. Next, the location of vessel 
information, HHH sub-band coefficients of the 3D discrete wavelet transformation (DWT), 
was detected by the proposed vessel-texture discrimination algorithm. Based on the initial 
segmentation, 3D DWT integrated with the 3D neutrosophic transformation could accurately 
detect the coronary arteries. 

Results: In order to evaluate the proposed method, the CTA scans of 20 patients were 
segmented for coronary arteries. The slice thickness was 0.625 mm, and the size of the 
volume was 512*512*(- ). Each sub-branch of the segmented coronary arteries was 
segmented correctly. To evaluate the performance of our method in segmenting coronary 
arteries, we compared our results with those ground truth values obtained from the 
commercial software from GE Healthcare and the level-set method proposed by Yang et al. 
[11]. We used the overlapping metric (OM) and Hausdorff distance ( d u ) to analyze the 
efficiency of each method and the results showed that the proposed method was much better. 

Conclusion: This work proposed a novel segmentation method for coronary arteries in CTA 
data sets. Based on the initial segmentation of coronary arteries obtained from 3D region 
growing, one-level 3D DWT and 3D neutrosophic transformation can be applied to detect 
coronary arteries accurately. In our experiment, we tested 20 data sets with delineated ground 
truth centerline trees. Results indicate that the proposed method is effective and useful It 
allows the coronary pathologies to be detected automatically and accurately. 

Index Terms — coronary arteries, 3D region growing, 3D discrete wavelet transformation, 3D 
neutrosophic transformation 



2 




I. Introduction 



Efficient and automatic image segmentation methods are useful for the isolation and 
visualization of vessels in computed tomographic angiography (CTA). There are many 
proposed method for the segmentation of vessels [1-14]. A vessel filter [1] can be used to 
enhance tubular structure, however, it cannot address the problem of the image force and 
veins, which can lead to a narrowed or broken segmentation of vessels. Parametric shape 
models [2-5] do not directly allow for the detection of topological changes, and they usually 
obtain a seriously narrowed segmentation in the neighborhood of a branch point in the vessel. 
Level-set approaches [6-13] are computationally expensive. They also suffer from leakage at 
places where the intensity gradients of the edges are relatively weak, and are very sensitive to 
the placement of the initial contour of the propagating front. Metz et al. [14] used the 
minimum cost path of the specified start and end points in vessel to detect the coronary 
arteries centerline. This is not an automatic method, detecting branches is difficult. Friman 
[15] proposed multiple hypothesis template tracking, which follows the direction of centerline 
obtained in advance. However, it is difficult to detect small branches and vessels by using this 
method. 

In this study, we propose a new method for automatically and correctly segment coronary 
arteries from CTA data sets. In image preprocessing, we detected the aorta automatically by 
using methods proposed in literature [5, 16]. The proposed coronary arteries segmentation 
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method is summarized as follows. First of all, we automatically obtained the seed point of a 



3D region growing by the difference between the two adjacent slices due to the small changes 
of the aorta between two adjacent slices. Next, 3D region growing was applied to initially 
search for the probable location of coronary arteries, which was then dilated by 3 voxels. 
Based on the dilation of the probable location, we detected the coronary arteries accurately by 
applying the 3D discrete wavelet transformation (DWT) and 3D neutrosophic transformation 
to the CTA volume. The location of vessel information, in HHH sub-band coefficients, was 
detected by the proposed vessel-texture discrimination algorithm. Accordingly, HHH 
sub-band coefficients were used, which were characterized and classified by a-means 
operation and K-means clustering. Finally, the proposed method was tested on several CTA 
data sets, and the experimental results indicated that the proposed method had a good 
performance. 

The rest of this study is organized as follows. Section II reviews some preliminaries, and 
section III uses 3D region growing and 3D DWT to propose a new method for segmenting 
coronary arteries. Section IV contains the experiments and discussion, and the conclusions are 
drawn in section V. 



II. Preliminaries 

In this section, we will briefly introduce the concepts of DWT and provide an overview of 
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some fundamental mathematical concepts that are used in this study. 



A. Region Growing 

Region growing is a simple, well-developed, region-based image segmentation 
technique [17]. It postulates that neighboring voxels within the same region have similar 
intensity values, and is also classified as a voxel-based image segmentation method since it 
involves the selection of initial seed points. In other words, this method of segmentation 
exa min es neighboring voxels of initial seed points and deter min es whether neighboring 
voxels should be added to the region. Consequentially, the general concept of region growing 
is to group voxels with the same or similar intensities to one region according to the given 
seed points and a homogeneity criterion. 

B. Discrete Wavelet Transform 

Wavelet transform is obtained by a single prototype function y/(x) which is regulated 
with scaling and shift parameters. To constmct y/{x ) , a scaling function <p(x) is determined. 
The discrete normalized scaling and wavelet basis functions are defined as 



<Pi,X 0 = : 


Pep ( t-n 


( 1 ) 


WiJ 1 )= 


(t-n 


( 2 ) 



where i and n are the dilation and translation parameters. Orthogonal wavelet basis 
functions not only provide a simple method to calculate coefficient expansion, but also span 
L 2 (\ i ) in signal processing. As a result, signal S(t ) e L 2 (l ) can be expressed as a series 
expansion of orthogonal scaling functions and wavelets. More specifically, 
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coefficients, respectively; j 0 is an integer to define an interval on which S(t) is a 
piecewise constant. The two-scale equations for scaling and wavelet basis function are given 
as follows. 

(Pit) = 72 Yj h ,„(Pi 2t ~ m ) 

me □ 

wit) = JiY.gJpilt-m) ( 5 ) 

meD 

where g m = (~\)" l h l . The coefficient h m in equation (4) has to meet several conditions for 

the set of the wavelet basis function to be unique, ortho normal, and have a certain degree of 
regularity. 

The coefficients h m and g m play a very crucial role in a given DWT. Performing the 
wavelet transformation does not require the explicit forms of (p(t) and y/(t) but only 
depends on h m andg m . The final output of the wavelet decomposition includes a set of 
j-level wavelet coefficients. One method to implement DWT is to use a filter bank that 
provides perfect reconstruction. DWT involves local analysis of frequency in space and time 
domains, and it provides multi-scale image details step by step. If the scale becomes smaller, 
every part becomes more accurate, and ultimately all imaging details can be focalized 
accurately. If DWT is applied to a volume, it will produce highest-frequency, 
middle-frequency, and lowest-frequency parts. Figure 1 shows the results of applying 3D 
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DWT to a volume, which includes eight parts: LLL, LLH, LHL, LHH, HLL, HLH, HHL, and 
HHH. The lowest-fiequency and highest-frequency parts are LLL andHHH respectively [16, 
17]. 
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Figure 1. The structure of applying 3D DWT to a volume. 



C. K-means Clustering 



K-means clustering is a method of cluster analysis which aims to partition n observations 
into k clusters in which each observation belongs to the cluster with the nearest mean. Given a 
set of observations (x l ,x 2 ,---,x n ), where each observation is a d-dimensional real vector, 
&-means clustering aims to partition the n observations into k ( k<n ) sets {iv, , .s 2 , - - \ so as 

to minimize the within-c luster sum of squares 



k 

arg min^ £ || x } 

i = 1 Xj GSj 




(6) 



7 












where //. is the mean of points in s f . 



III. The Proposed Segmentation Method 

In order to segment coronary arteries accurately from CTA data sets, 3D region growing 
was initially applied to search for the probable location of the coronary arteries. Next, we 
used the 3D DWT and 3D neutrosophic transformation to accurately detect the coronary 
arteries. 

A. Initial Segmentation of Coronary Arteries 

This section discusses the initial segmentation of coronary arteries using 3D region 
growing. In image preprocessing, we found the aorta automatically by using methods 
proposed in literature [5, 16]. The selection of the seed point of 3D region growing was 
initially made to check which slice began the information of the coronary arteries. Due to the 
small changes in aorta area between the two adjacent slices, we automatically obtained the 
seed points by the difference between the two adjacent slices. 

As shown in Figure 2, we use the difference between the two adjacent slices (a) and (b) to 
automatically obtain the seed points bounded by the blue line in (c) which indicates the 
boundary of coronary arteries. The boundary of coronary arteries denotes the high-frequency 
sub-band d j fk) in equation (3) when comparing vessel lumen and background. 





(C) 

Figure 2. The 3D region growing seed points in (c) were automatically obtained by using the 
difference between the two areas bounded by the blue lines in adjacent slices (a) and (b). 

Since coronary arteries do not exhibit abrupt intensity changes along their centerline [4], 
a rough tubular mask of coronary arteries can be easily constructed by 3D region growing. 
We chose a 26-connected neighborhood for our adjacent pixel relationship, and then the 3D 
region growing method was applied with a set of pre -specified seed voxel(s) and grown from 
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these seeds by merging neighboring voxels whose properties were most similar to the 



pre-merged region. The homogeneity criterion was defined as the difference between the 
intensity of the candidate voxel and the average intensity of the pre-merged region. The 
selection of the seed point was initially intended to check which slice began the information 
of the coronary arteries. Next, the homogeneity criterion was applied to group voxels with the 
same or similar intensities into one region. If the homogeneity criterion was satisfied, the 
candidate voxel was merged with the pre-merged region. The process was repeated until no 
more voxels were assigned to the region, and then the number of all merged voxels was 
calculated. In order to avoid leakage, the total number of merged voxels was limited to 12000 
otherwise, 3D region growing was restarted by automatically using an improved homogeneity 
criterion. Finally, the initial segmentation of the coronary arteries in a volume was completed. 

B. Vessel-Texture Discrimination 

According to Parseval’s theorem, the energy in a signal S(t) is given as follows [16, 17] : 

co oo co 2 

j \s(t )\ 2 dt=x |c (/)| 2 + z Z \ d M <?> 

/=- 00 j = 0 k=— 00 

This equation implies that the energy of a signal is the summation of low-frequency and 
high-frequency coefficients. DWT is a good analytic tool for image texture analysis or 
line -based patterns [21-24], Since a vessel is a type of 3D line -based pattern in CT volume, 
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we used l 2 -norm of these wavelet coefficients to find the energy of line-based patterns. 



l 2 -norm was defined as 

£( cH | c || 2 = 2>,.| 2 ( 8 ) 

/=i 

where the vector C = [c i ] lx was the wavelet coefficients of a frequency channel. The 
searching algorithm is summarized as follows: 

Algorithm 1 : 

1. Transform a given vessel volume into frequency channels by a specified number of 
decomposition levels. We usually set the number to one in the first search. 

2. Use equation (8) to calculate the averaged / 2 -norm of each channel and maximum of 
E(C) lor the vessel volume. 

3. If the maximum of E(C) was significantly greater than another channel’s E(C ) , the 
search was stopped. Otherwise, the number of decomposition levels was increased 
followed by a repeat of stepl. 

By using the above algorithm, we observed that the most significant information of the 
vessel texture often appeared in the high frequency channels. Thus, we used the sub-band 
HHH to detect vessels in this study. 

C. Accurate Detection of Coronary Arteries 
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The initial segmentation of coronary arteries in a volume was completed by 3D region 



growing as described in sub-section A. Since region growing is a simple region-based image 
segmentation method, it was only used to search for the initial location of the coronary 
arteries. We then accurately detected the coronary arteries by applying DWT to each slice in 
the volume, as described in this sub-section. 



First, the initial location of the coronary arteries obtained from 3D region growing was 
dilated by 3 voxels. Next, we used Haar wavelet bases in (1) and (2) to transform the host 
images into the orthogonal DWT domain by one -level decomposition. Only HHH sub-bands 
were employed for further processes, because most of the information on the coronary arteries 
and boundaries were in the HHH sub-bands. We calculated the mean energy using 
coefficients of HHH sub-bands in a local window w as follows: 



where 



1 r s 1 

HHH (i, j, k) = ZEE HHH (l, m , n ), 

WXWXW l=i m=jn=k 



(9) 
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Next, the sub-bands HHH was characterized by 3 membership sets T , F , and U . 
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where 



HHH mi „ = min {////// (i, ./,/:)} , 



HHH mm =max.{HHH(i,j,k)}, 



S(i, j,k) = \HHH(i, j, k ) - HHH (i, j,k) I, 



f>„„n = niin ./U' ) }, 



$■* = max {<?(/, 7,*)}. 



That is, a pixel could be represented as a neutrosophic domain P{t,f,u) which means the 
pixel is t% true, f% false, and u% uncertain, where t varies in T , f varies in F , 
and u varies in U ■ In order to reduce the uncertainty u% , a-means operation was 
employed as follows. 
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The inverse DWT was then applied to obtain a new volume which possessed the true subset. 
Finally, we applied K-means clustering (K=3) in (3) to differentiate vessel lumen, vessel 
boundary (true subset) and background. The true subset T was retained respectively. 



IV. Experiments and Discussion 



To test the proposed method, CTA volumes obtained from a CT system were segmented 
for coronary arteries. The slice thickness was 0.625 mm and the volume was 512*5 12*( • ) in 
different data sets. The window size w was set to 3 which was enough to capture the local 
texture characteristics. The parameter a was set to 0.2. We tested 20 data sets, most of 
which were segmented successfully except that for a few small branches that were lost in 2 of 
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the data sets due to the local failure in region growing. To evaluate the performance of our 
segmented coronary arteries, we compared our results with that obtained from the ground 
truth values obtained from the commercial software from GE Healthcare and the level-set 
method. We used the overlapping metric (OM) and Hausdorff distance (d n ) to analyze the 
efficiency of each method. 

A. The Segmenting Efficiency of 2D Imaging 

Figure 3(b) shows the results of 6 slices in 1 CTA volume obtained using the proposed 
method. The areas bounded by the red line are coronary arteries. The segmenting efficiency 
was compared with the manually delineated ground truth data n r in Figure 3 (a) by using 
OM which was defined as 

JN T nNA 

OM =2 -T ^ 

V N T + , 

where n t indicates the pixels/voxels of the segmented coronary arteries. The OM was close 
to 1 when the segmentation was well matched to the reference ground truth, and approached 
zero when the results had no similarity to the reference. 

In the 6 slices in Figure 3(b), the segmentation results showed that the proposed method 
detected coronary arteries accurately. As shown in Table I, the average OM of the proposed 
method was 0.96. 
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(a) ground truth data (b) the segmented results 



Figure 3. Comparison of the 2D segmentation (b) with respect to the ground truth data (a). 

TABLE I 

Comparison of the 2D Image Segmentation Results. 





Yang etal. [11] 
Method 


Proposed 

Method 


Mean OM 


0.68 


0.96 



B. The Segmenting Efficiency on a 3D Volume 



The first focus of the comparison was the correctness of the 4 main branches: the right 
coronary artery (RCA), the left anterior descending artery (LAD), the circumflex (CRX), and 
the first diagonal artery (DA). Figures 4 and 5 show the detected coronary arteries in the CTA 
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volume obtained from the GE Healthcare and the proposed method, and the proposed method 



segmented these 4 main branches correctly compared to the 4 main branches in Figure 5. 

Another test of performance is the correctness of the remaining branches. Due to the 
multi-resolution of the DWT, each sub-branch in Figures 3 and 5 was correctly segmented. 
Furthermore, the coronary arteries obtained from the proposed method were much better than 
those obtained from the level-set method, which had several leakages as seen in Figure 6. 

In addition to the OM, the difference between the segmented vessel surface and the 
manually delineated ground truth data was measured by the Hausdorff distance [25] which 
was defined as follows: 

cl H {X,Y ) = maxi sup inf d(x, y), sup inf d(x,y ) 

[xeX y sY ysY AeX 

where X and V are the vertices of the mesh surfaces of the arteries corresponding to the 
segmentation results and the ground truth and d(x,y) measure the Euclidean distance 
between points x and y belonging to vertices X and Y . Table II lists the mean OM and 
mean Hausdorff distance for the proposed method and the method of Yang et al. [11]. The 
results show that the proposed method was much better than that of Y ang et al. [11] in terms 
of both OM and Hausdorff distance. 
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C. Diameter Measurement 



In this sub-section, we computed the diameters of segmented coronary arteries using the 
proposed method and the method of Yang et al. Many efficient algorithms have been 
proposed to extract of the tube centerline. We applied the algorithm proposed in by Lee et al. 
[26] to extract the centerline of the segmented coronary arteries. Using these extracted 
centerlines, we obtained the cross sections of the segmented coronary arteries, as shown in 
Figure 7. By computing the area A of each cross section, the diameter r was estimated as 
follows: 



r 

Table III shows the estimated diameter of the same cross section for the proposed method 
and the method of Yang et al. [11]. The diameter of the proposed method was closer to the 
diameter in ground truth data at the same cross section than that obtained from the method 
Yang etal. 
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Figure 4. The 3D coronary arteries manually edited by a radiologist using an AW workstation 



(GE Healthcare, Wisconsin, USA). 




Figure 5. The 3D coronary arteries in the CTA volume obtained from the proposed method. 
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Figure 6. The 3D coronary arteries in the CTA volume obtained from the level set approach. 




(a) 



21 




(b) 



Figure 7. (a) The cross section of the segmented coronary arteries, (b) Magnified view of 



(a). 



TABLE II 

Comparison of the 3D Segmentation Results. 





Yang etal. [11] 
Method 


Proposed 

Method 


Mean OM 


0.60 


0.92 


Mean d H (X,Y ) 


1.08 


0.68 



TABLE III 

Comparison of the Cross-sectional Diameter. 



Ground Truth 


Yang etal. [11] 
Method 


Proposed 

Method 


4.7 


4.4 


4.8 
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D. Experimental Environment and Execution Time 



The proposed method was implemented in MATLAB (R2011a) on a standard 
specification PC with a 3.2 GHz CPU and 12 GB RAM. The average execution time was 58 
seconds to extract the entire coronary tree, compared to approximately 47 seconds for the 
method by Yang et al. for the same process. 

V. Conclusions 

Accurate extraction of coronary arteries is important to assess artery lesions in clinical 
practice. In this study, we propose a novel method to segment coronary arteries automatically. 
Based on the initial segmentation obtained from 3D region growing, one-level 3D DWT and 
3D neutrosophic transformation were applied to detect coronary arteries accurately. The 
location of vessel information, in HHH sub-band coefficients of DWT, was successfully 
detected by the proposed vessel-texture discrimination algorithm. Accordingly, the HHH 
sub-band coefficients were used and characterized and classified by 3D neutrosophic 
transformation and K-means clustering. The experimental results verify the efficiency of the 
proposed method. 
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